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Abstract 


We explore a response of a non-linear non-axisymmetric mean-field solar dy¬ 
namo model to shallow non-axisymmetric perturbations. After a relaxation pe¬ 
riod the amplitude of the non-axisymmetric held depends on the initial condition, 
helicity conservation, and the depth of perturbation, ft is found that a pertur¬ 
bation which is anchored at O.9i? 0 has a profound effect on the dynamo pro¬ 
cess, producing a transient magnetic cycle of the axisynnnetric magnetic held, 
if it is initiated at the growing phase of the cycle. The non-symmetric with 
respect to the equator perturbation results in a hemispheric asymmetry of the 
magnetic activity. The evolution of the axisynnnetric and non-axisymmetric 
helds depend on the turbulent magnetic Reynolds number R m . In the range of 
R m = 10 4 — 10 6 the evolution returns to the normal course in the next cycle, in 
which the non-axisymmetric held is generated due to a non-linear a-effect and 
magnetic buoyancy. In the stationary state the large-scale magnetic held demon¬ 
strates a phenomenon of “active longitudes” with cyclic 180° “hip-hop” changes 
of the large-scale magnetic held orientation. The hip-hop effect is known from 
observations of solar and stellar magnetic cycles. However, this effect disappears 
in the model which includes the meridional circulation pattern determined by 
helioseismology. The rotation rate of the non-axisymmetric held components 
varies during the relaxation period, and carries important information about the 
dynamo process. 
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1. Introduction 


Dynamo theories commonly assume that the magnetic activity of the Sun is 
approximately axisymmetric on large spatial (size of the Sun) and temporal (the period of 
solar cycle) scales. These models provide a quantitative self-consistent description of the 
22-year solar magnetic cycles, and allow us to investigate the basic mechanisms of the solar 
dynamo. However, deviations from the axisymmetry are rather strong at any particular 
moment of observations. Intermittent patterns of magnetic fields on the solar surface are 
formed because the magnetic held emerges on the surface like separated magnetic patches, 
e.g, in the form of sunspot groups. Such phenomena make a significant contribution to the 
large-scale non-axisymmetric magnetic held of the Sun. 


Raedler (1986) discussed dynamo generation of large-scale non-axisymmetric (NA) 


magnetic held on the Sun. It was found that the differential rotation suppresses generation of 
the NA magnetic held. Non-linear dynamo processes can maintain a weak non-axisymmetric 


held in expense of the axisymmetric (AS) magnetic held (see, Raedler et al. 1990, Moss 


1999, Elstner & Korhonen 2005). Also it was found that the dynamo generated NA 


magnetic held rotates rigidly. Further theoretical developments (see, e.g., Moss et al. 


1991, Bigazzi & Ruzmaikin 2004; Moss 2004, Berdyugina et al. 2006) showed that some 


of the properties of the non-linear non-axisymmetric mean-held dynamo can be invoked 
for interpretation of the origin and evolution of the so-called “active-longitudes” of solar 
magnetic activity. Phenomenon of the active longitudes (AL) is probably one of the most 
interesting manifestations of the solar non-axisymmetric magnetic held. It appears when 
the solar activity persists within the hxed interval of longitudes for a long period of time 


(see e.g., Vitinskii 1966, Burnba & Howard 1969, Vitinsky et al. 1986). 


We have to mention that the question about persistence of the AL on century 
time scales interval remains a highly controversial issue both from the observational and 
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theoretical points of view. For example, Berdyugina (2004), Berdyugina et al. (2006) and 


Zhang et al. (2011) reported about the AL which are persistent over a century long time 


interval. However, Pelt et al. (2010) found that the AL have the maximal lifetime of about 


one solar cycle. Another phenomenon which is related to the AL is the so-called “hot spots” 


of the solar flare activity (Bai 1987, 2003). The reason for the different name is because 


in general the longitudinal position of the activity nests is different in the Northern and 
Southern hemispheres. Bai (2003) found that those “hot spot could persist rigidly rotating 
with period about 27 days up to three solar cycles. The nonlinear dynamo models predict 
that the energy of the non-axisymmetric modes is only about 10 -4 of the energy of the 
axisymmetric (AS) component (Berdyugina et al. 2006). It is not clear how such a weak 
magnetic field can modulate the nests of the sunspots activity. 

While, the standard mean-field model can not explain the origin of the NA magnetic 


field (Raedler 1986), it was noted that the solar activity produces the large-scale NA 


modes which are well seen in the coronal hole configurations (Stix 1977). The very recent 


example of such events was observed by the SDO/HMI during the last decade of May 


2015. Observations at the Wilcox Solar Observatory (Duvall et al. 1979, Hoeksema 1995) 


found that the strength of the NA modes of the radial magnetic field, e.g., the mode with 
the azimuthal number m=l, can be about 1 G during epoch of the solar maxima. The 
axisymmetric dipole has the same magnitude during the solar minims. It is likely that the 
origin of this NA held is related to decay of solar active regions. However it is unclear how 
the evolution of such NA held may impact the solar dynamo process. The effect of the NA 
held on the global dynamo has not been studied before. 


In this paper we explore a non-linear response of a mean-held dynamo model to a 
shallow NA m=l perturbations with the held strength of 1G. It is assumed that these 
perturbations result from abrupt instability-type events, which can be described by injecting 
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NA perturbations into the system over a very short period of time. We consider a fairly 
complete theoretical description of the mean turbulent electro-motive force taking into 
account the known properties of the solar convection zone and including the anisotropic 
turbulent effects due to the global rotation. The model includes a nonlinear magnetic 
buoyancy effect, and two types of non-linearity in the a-effect, described as “algebraic” and 


“dynamical” quenching (Tworkowski et al. 1998). The algebraic quenching is due to the 
back-reaction of the dynamo-generated magnetic held on helical turbulence. The dynamical 


quenching results from a magnetic helicity conservation condition (Kleeorin & Ruzmaikin 


1982). We will show that the nonlinear non-axisymmetric effects are sufficiently strong to 


reproduce the “hip-hop” phenomenon and explain the rotation rate of active longitudes on 


the Sun (Tuominen et al. 2002a, Berdyugina et al. 2006, Gyenge et al. 2012). 


2. Basic equations 

Evolution of the large-scale magnetic held in perfectly conductive media is described 


by the mean-held induction equation (Krause & Radler 1980; Moffatt 1978, Parker 1979): 


d t (B) = V x (£+ (U) x (B» 


( 1 ) 


where £ = (u x b) is the mean electromotive force; u and b are the turbulent fluctuating 
velocity and magnetic held respectively; and (U) and (B) are the mean velocity and 
magnetic held. Our solution of the dynamo equation will follow the outline given earlier by 


Moss et al. (1991) and Moss (1999). For convenience we decompose the magnetic held into 


the axisymmetric (AS), (hereafter B-held), and non-axisymmetric (NA) parts, (hereafter 
B-held): (B) = B + B. We assume that the mean how is axisymmetric (U) = U. Let 
(p = and r = re r be vectors in the azimuthal and radial directions respectively, then we 
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represent the mean magnetic field vectors as follows: 


(B) 

= B + B 

(2) 

B 

= (j)B + V x 

(3) 

B 

= V x (f T) + V x V x (f 5), 

(4) 


where A, B , T and S are scalar functions representing the AS and NA parts respectively. 
Assuming that A and B do not depend on longitude, Eqs([3j [4]) ensure that the field (B) is 
divergence-free. Taking the scalar product of Eq([lJ with vector <p we get equations for the 
AS magnetic field components, 


d t B = 0 • V x (£+U x B) , 
d t A = 4> ■ (£+U x B) , 


(5) 

( 6 ) 


To get equation for T we take curl of Eq([TJ) , and then calculate the scalar product with 
vector r. Similarly, equation for S is obtained by taking twice curl of Eq(JTj) and then the 
scalar product with vector r. The procedure is described in detail by |Krause fe Radler 


(1980). Equations for the NA field are 


d t A n T = AnK^ + AnK^, (7) 

d t A n S = A n uW + A n uW, (8) 


<9/i d^i sin 2 6 dtp 2 ’ 

A n V {u) = -r.VxVx(UxB), (9) 

A n V {£) = -r-VxVxS, (10) 

A n U {u) = -f ■ V x (u x B) , (11) 

A n U^ £) = -i'Vxf. (12) 


where A 
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The scalar functions with superscript (u ' > contain contributions from the large-scale AS flows 
like the differential rotation or meridional circulation. The integration domain includes the 
solar convection zone from 0.71 to 0.99 R®. The distribution of the mean flows is given by 


helioseismology (Howe et al. 2011 and Zhao et al. 2013). Profiles of the angular velocity 


and meridional circulation are illustrated in Figure 1. 

We use formulation for the mean electromotive force obtained by Pipin(2008). 

The calculations of the mean electromotive force are done using the mean-field 
magnetohydrodynamics framework and the so-called “minimal tau-approximation” (see, 
e.g., Blackman & Field 2002; Radler et al. 2003; Brandenburg & Subramanian 2005[ ). The 
tau-approximation suggests that the second-order correlations do not vary significantly on 
timescale r c which corresponds to a typical turnover time of the convective flows. The 
theoretical calculations are performed for anelastic turbulent flows. They take into account 
the effects of density stratification, spatial inhomogeneity of the intensity of turbulent 
flows, and inhomogeneity of the large-scale magnetic fields. The effects of large-scale 
inhomogeneity of the turbulent flows and magnetic fields are computed in the first order 
of the Taylor expansion in terms of ratio £/L, where £ is a typical spatial scale of the 
turbulence, and L is a spatial scale of the mean quantities. The mean electromotive force, 
£, is expressed as follows (Pipin, 2008): 


£i — {oiij + 7 ij) {B)j TJijkVj (B) k . (13) 

where symmetric tensor a t j models the generation of magnetic field by the a- effect; 
antisymmetric tensor 7 ^ controls the mean drift of the large-scale magnetic fields in 
turbulent medium; tensor rj t jk governs the turbulent diffusion. We take into account the 
effect of rotation and magnetic field on the mean-electromotive force (see, e.g., 


Pipin 2015 


for details). To determine unique solution of Eqs.([5]|8]) we apply the following gauge (see, 




























e.g., Krause & Radler 1980, Bigazzi & Ruzmaikin 2004): 


-2n r l 


»2-k r 1 


Sdjidp = 0, 


Td/id(j> = 0 . 


'o J -1 


'0 J -1 


The same gauge is used in Eqs(|9 12). 


(14) 


2.1. Nonlinear interaction of the axisymmetric and non-axisymmetric modes 


Interaction between the AS and NA modes in the mean-field dynamo models can be 


due to nonlinear dynamo effects, for example, the o-effect, (Krause & Radler 1980, Moss 


1999). In our model the a effect takes into account the kinetic and magnetic helicities in 


the following form: 


a 




= C a sin 2 9i/j a (/?) qA j t]t + aA 


(H) (M) (X') T c 


(15) 


lJ 4n pP 

where C a is a free parameter which controls the strength of the a- effect due to turbulent 
kinetic helicity; and express the kinetic and magnetic helicity parts of the 
a-effect, respectively; tjt is the magnetic diffusion coefficient, and (y) = (a • b) (a and 
b are the fluctuating parts of magnetic field vector-potential and magnetic field vector). 
Both the a-A and depend on the Coriolis number fl* = 47T—4—, where P rot is the 

* rot 

rotational period, r c is the convective turnover time, and i is a typical length of the 
convective flows (the mixing length). A theoretical justification for the latitudinal factor, 
sin 2 9, in Eq(15) was given by Kleeorin & Rogachevskii (2003). Function i/j a (p) controls the 


so-called “algebraic” quenching of the a- effect where (3 = | (B) | / \J4irpu' 2 , u' is the RMS of 
the convective velocity. For the case of the strong magnetic field, (3 1^*1, ip a P~ 2 . The 
interaction between the axisymmetric and non-axisymmetric dynamo modes via ?/; Q (/9) is 
because both modes contribute to parameter {3. Also, for the case [3 >1, the latitudinal 
profile of the a effect changes. This can affect the dynamo conditions for excitation of the 
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NA modes. Raedler et al. (1990) and Moss (1999) discussed the evolution of NA magnetic 


field in a simple dynamo model with such “algebraic” quenching. 


The dynamical quenching is caused by the magnetic helicity conservation (see, Kleeorin 


& Ruzmaikin 

1982 

). This effect was discovered by 

Frisch et al. 

(1975 

) and Pouquet et al. 


(1975). Contribution of the magnetic helicity to the a-effect is expressed by the second 


term in Eq.(15). The magnetic helicity density of turbulent field, (y) = (a • b), is governed 


by the conservation law (Hubbard & Brandenburg 2012, Pipin et al. 2013): 


d(x) 


(tot) 


(x) 


-277(B)-(J)-V.^, 


(16) 


dt R m r c 

where (y)^°^ = (x) + (A) • (B) is the total magnetic helicity density of the mean and 
turbulent fields, 7- x = — (y) is the diffusive flux of the turbulent magnetic helicity, and 

R m is the magnetic Reynolds number. The coefficient of the turbulent helicity diffusivity, 


r) x , is chosen ten times smaller than the isotropic part of the magnetic diffusivity (Mitra 


et al. 2010): rj x = Similarly to the magnetic field, the mean magnetic helicity 


density can be formally decomposed into the axisymmetric and non-axisymmetric parts: 
(x) < ' tot ' 1 = yb°0 -f yh°d. The game can be done for the magnetic helicity density of the 


turbulent field: (y) = y + y, where y = a ■ b and y = (a ■ b). Then we have, 


y( tot ) = y + A-B + A-B, 


X 


(tot) _ 


= y + AB + AB + AB, 


(17) 

(18) 


Evolution of the y and y is governed by the corresponding parts of Eq(Jl6|). Thus, the model 
takes into account contributions of the AS and NA fields in the whole magnetic helicity 
density balance, providing a non-linear coupling. We see that the a-effect is dynamically 
linked to the longitudinally averaged magnetic helicity of the NA B-field, which is the last 


term in Eq(17). Thus, the nonlinear a-effect is non-axisymmetric, and it results in coupling 


between the AS and NA modes. The coupling works in both directions. For instance, the 
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azimuthal a-effect results in (B^). If we denote the NA part of the by 

then the mean electromotive force is £^ + a^B^. This introduces a new source in 

Eq(| 6 ]) which is usually ignored in the axisymmetric dynamo models. 


Magnetic buoyancy is another nonlinear effect which is important in the large-scale 
dynamo. The part of the mean electro-motive force which is responsible for magnetic 


buoyancy is (Kichatinov & Pipin 1993): 


£ {i3) = Vgr x B, 


(19) 


where Vg = Cg MLT — /3 2 K (j3 ), u! is the RMS convection velocity, K (j3) can be found 
7 

in (Kichatinov & Pipin 1993; Pipin 2015), 7 is the adiabatic exponent, olmlt is the 
mixing-length theory parameter, Cg is a free parameter to switch on/off this effect in the 
model. For the case /) C 1, ~ 1 and the up-flow velocity, Vg, is proportional to the 

pressure of large-scale magnetic field. Similarly to the a-effect, the Vg is non-axisymmetric 
and contributes to the source terms in Eqs( [5p| ). Note that advection of the large-scale 
magnetic field by the magnetic buoyancy reduces concentration of the magnetic field near 
the bottom of the convection zone, and increases the field strength near the top. 


2.2. Parameters of the convection zone and numerical procedure 


The distribution of the turbulent parameters, such as the typical convective turn-over 
time, t c , the mixing length, £, and the RMS convection velocity, u', are taken from the 
solar interior model of 


Stix 


(2002). We define the mixing-length: £ = Omlt I I ', where 


= V log p is the inverse pressure scale, and the mixing-length parameter Omlt = 2. 


The profile of the turbulent diffusivity is taken in the form t]t = C Tl 


u , 2 t c 


3 fa 


-, where 


f ov (r) = 1 + exp [50 (r ov — r)], r ov = 0.725 R Q controls quenching of the turbulent effects near 
the bottom of the convection zone, which is rg = 0.715A. Free parameter C v , (0 < C v < 1) 
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controls the efficiency of mixing of the large-scale magnetic field by turbulence. It is usually 
employed to tune the period of the dynamo cycle. 

The numerical scheme employs the spherical harmonics decomposition for the 
non-axisymmetric part of the problem, i.e., the scalar functions T and S in Eqs( [7p| ) 
represented in the form: 


are 


T (r, n,<j>, t) = N ij, m (r, t) P, H exp , 


5 (r, p, 0, t) = S L , m (r,t) p/ m| exp (im<p ), 


( 20 ) 

( 21 ) 


where P/ n is the normalized associated Legendre function of degree l > 1 and order m > 1. 
The simulations which we will discuss include 600 spherical harmonics [l m ax = 28). Note 
that S(_ m = Sf m and the same for T. We employ the pseudo-spectral approach for 
integration along latitude. The second-order finite differences are used for discretization int 
the radial direction. The numerical integration is carried out in latitude from the pole to 
pole and in radius from r& = O.715P 0 to r e = 0.99-Rq. All the nonlinear terms are calculated 
in the real space. The transformation between the spectral spherical harmonic and the 
real 3D space was done using the Intel Fortran FFT library. We implement algorithms 


of Muciaccia et al. (1997) to speed-up the transform calculations. At the bottom of the 


convection zone we set up a perfectly conducting boundary condition for the axisymmetric 
magnetic field, and for the non-axisymmetric field we set the functions S and T to zero. 
At the top of the convection zone the poloidal held is smoothly matched to the external 
potential held. The boundary conditions for toroidal held allow held penetrate the surface 


(Moss et al. (1991), Pipin & Kosovichev (2011)): 


5 r ^B + (l-8)S e = 0, 
- 


( 22 ) 


0 


(23) 
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where parameter 6 = 0.99. 


The particular choice of parameters was discussed in our previous papers (see, e.g, 


Pipin & Kosovichev 2014). The free parameters are C a = 0.04, Cs = \C a , C v = -^ and 


the anisotropy parameter a = 3 (see Pipin & Kosovichev 2014). The a-effect parameter 
C a is about 30% above the dynamo generation threshold. For the chosen values of C n and 
a, the turbulent diffusion coefficient in the near-surface shear layer, at r = O.9i? 0 is about 


10 9 m 2 s 1 which is in agreement with surface observations (see, Abramenko et al. 2011). 


The magnetic helicity conservation is determined by the magnetic Reynolds number R m , 
for which we considered values: 10 4 - 10 6 . 


To investigate the influence of meridional circulation we consider two models: Ml 
without meridional circulation and model M2 with the double-cell meridional circulation 
with a characteristic velocity 10 m/s. For the model M2 we employ the larger parameter C a 


because this model has a larger critical dynamo threshold (see, Pipin & Kosovichev 2013 


Pipin & Kosovichev 2014). For this value of C a model M2 has the same magnitude of the 


generated AS toroidal magnetic held inside the convection zone as the model Ml. The set 
of parameters in the models is summarized in Table 1. 


To quantify the mirror symmetry type of the toroidal magnetic held relative to the 
equator we introduce the parity index, P , as follows, 


P = (24) 

E s ,a = 1 f (B (r s , 0) ± B(r s , -0)) 2 sin0d0, 


where “+” corresponds to the index “S”, and ” is for “A”; E = Es + Ea, Es and Ea are 
energies of the symmetric and anti-symmetric components of the AS toroidal magnetic held 
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at r s = 0.9-Rq . Similarly we define parameter P for the NA magnetic field, 

Es — Ea 


P = 
Es,a = 


(25) 


B 0 (r s , 9 , </>) ± B ^ (r s , -9, </>)) sin Odddcj), 


where E = E s + Ea- The parity of the total magnetic held is P = ( PE + PE j / [E + Ej. 




We will see that perturbation can affect the cycle amplitude. Similar to Raedler et ah 


(1990) we introduce a parameter to measure deviation of the toroidal magnetic held from 


symmetry about the axis of rotation 


M — 1 — — 


E 


E + E 


(26) 


We simulate the sunspot number W using the anzatz (Pipin et al. 2012): 

W ( t) = (B max ) exp y ) > ( 27 ) 

where in model Ml (-B max ) is the maximum strength of the toroidal magnetic held averaged 
in the subsurface layers over radius in the range of 0.9 — 0.99 R & , and Bq is a characteristic 
strength of the toroidal magnetic held, Bq = 800G. In model M2 we measure (B max ) in 
the layer of convergence of the two meridional circulation cells, which is in the range of 
0.85 —O.9i? 0 . This layer corresponds to the maximum of the toroidal magnetic held strength 
in the convection zone for model M2. The AS model with the double-cell meridional 


circulation was discussed in details in our previous papers (Pipin & Kosovichev 2013, Pipin 


& Kosovichev 2014). 


2.3. Initial conditions 

In our hrst runs the weak initial held, which consisted of a superposition sum of polar 


and equatorial dipoles with the magnetic held strength of 0.01G, evolved to a state in 
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which the axisymmetric dynamo regime dominates. In this regime, the typical strength 
of the axisymmetric toroidal field in the convection zone is about lkG. However, the 
non-axisymmetric field is rather weak with the strength about 10 _5 G. This means that 
in our model the non-axisymmetric magnetic field is linearly stable, unless the dynamo 


governing parameters are forced to be non-axisymmetric, e.g., like in the paper by Bigazzi 


& Ruzmaikin (2004). 


Exploring the nonlinear solutions we found that the evolution of the non-axisymmetric 
field depends on the initial conditions which include the strength and geometry of both 
the axisymmetric B-field and non-axisymmetric B-field. The evolution of the large-scale 
magnetic field depends on the presence of the meridional circulation, too. 


In the following section we present results for the non-axisymmetric dynamo which was 
perturbed by a finite-amplitude non-axisymmetric B-field in the developed axisymmetric 
dynamo regime. Such non-axisymmetric perturbations can be developed either due to 
evolution of active regions or due to instabilities not described by the mean-field theory 


(e.g., Dikpati & Gilman 2001). For the seed field we consider a lion-symmetric relative to 
the equator perturbation represented by a sum of the equatorial dipole (1=1, m=±l) and 


quadrupole (/ = 2, m = ±1) components. In Eq(21) we define 


Si, i = — ( 1 — erf 


To 


Si ,2 = - ( 1 - erf 


d 

!r s - r) 


d 


r \ 2 

1 e 


(28) 

(29) 


where, r e = 0.99-Rq, r s = 0.9 R & is the bottom of the subsurface shear layer, d = 0.02 R®. 
The other Si, m and T) m coefficients are zero in the perturbation. The initial non- 
axisymmetric perturbation is concentrated in the near-surface shear layer. The depth of the 
non-axisymmetric perturbation can influence the evolution of the axisymmetric dynamo. At 
the moment of the initialization of perturbation the strength of the axisymmetric toroidal 
field is by two orders of magnitude greater than of the non-axisymmetric one. At the same 












15 


time the maximum amplitude of the NA radial magnetic field at the surface is about twice 
of the axisymmetric one (see Fig. 2c). Also, the initialization time we have the following 
parameters for the parity of the AS and NA parts are P = —1, P = —0.3 respectively and 
the indexes of the axisymmetry are M 10 -4 . 


3. Results 


Figure [2] illustrates the geometry of the axisymmetric and non-axisymmetric fields just 
before initialization of the non-axisymmetric perturbation. The perturbation is initialized 
at t ~ 13.2 years for model Ml and t ~ 5.7 for model M2. Those times corresponds 
approximately the same phase of the AS toroidal magnetic field in the upper half of the 
convection zone in the models. Figure [3] illustrates the evolution of the axisymmetric 
magnetic field before and after the perturbation for two models, with and without the 
meridional circulation. We show the time-latitude diagrams for the toroidal magnetic held 
in the subsurface shear layer and the radial magnetic held at the surface. In the model 
with the meridional circulation the toroidal magnetic held is shown for the middle of the 
convection zone (see 


Pipin & Kosovichev 2013). The radial evolution is shown for 30° 


latitude in the Northern hemisphere. 


The models show that the imposed non-axisymmetric perturbation produces a transient 
cycle in both models Ml and in M2. In the Northern hemisphere, where the initiated 
perturbation is greater, the simulated sunspot number cycle is stronger than in the Southern 
hemisphere. The perturbation affects the reversal of the polar magnetic fields. In the 
Northern hemisphere the polar field reversal occurs earlier than in the Southern hemisphere 
where we see multiple reversals. Model Ml shows a time shift of about 2 years for the polar 
reversals in the Northern and Southern hemispheres. We find that these phenomena depend 


on the depth of perturbation (parameters, r s and d in Eq( 28)). For instance, the polar 
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Fig. 1.— a) The isolines of constant angular velocity ranging from 0.6fio to 0.96h2o (flo = 
2.87 x 1CT 6 s _1 ); b) illustration of the double-cell meridional circulation model, consistent 


with the helioseismology results of Zhao et al. (2013). 


Table 1: Summary of the dynamo models and their parameters. 


Components 

Ml 

M2 

The S 

Eq. (13) 

same 

circulation 

no 

U = lOm/s 

free parameters, (see. 

Pipin 

2015 

) 

= 0.04, Cs = c a / 3, 

C v = 0.06, a=3, Cp = 1 

and Cp = C v , R m = 10 4 ~ 6 , 

OiMLT = 2 

C a = 0.05, others are same 

as in Ml 
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Fig. 2.— Axisymmetric and non-axisymmetric field structure at the moment of initialization 
of the non-axisymmetric perturbation: a) distribution of the axisymmetric toroidal magnetic 
held (background image) and poloidal held lines in the meridional cross-section for model 
Ml; b) the same as a) for model M2; c) components of the non-axisymmetric magnetic held 
at the surface, B r is shown as color image, B$ is shown by contour lines plotted every 3G; 
d) illustration of the magnetic held lines in model Ml at the moment of initialization of the 
perturbation. 


















18 



TIME, [YR] 



5 10 15 20 25 30 

TIME, [YR] 


Fig. 3.— . a-c) Dynamo evolution in model Ml before and after the initialization of the non- 
axisymmetric perturbation at t = 13 yr: a) time-latitude diagram, b) time-radius diagram at 
30° latitude, c) the simulated sunspot number for the Northern and Southern hemispheres; d- 
f) the same as in a-c) for model M2. Color images show the radial magnetic held component. 
The contour lines show the toroidal component, plotted every 100 G. 
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field reversal happens earlier in the Southern hemisphere than in the Northern one if the 
imposed perturbation is shallower r s = 0.95 R & . We also see that in both models the cycle 
returns quickly to the previously established axisymmetric state. Additional runs which are 
not illustrated here show that this restoration could take a longer time interval (more than 
1 cycle) if the axisymmetric field would have a mixed parity at the initialization moment. 

The evolution of the magnetic helicity density (Eq [l6| depends on the magnetic 
Reynods number, R m , which is a free parameter of the model. For higher R rn the relaxation 
in both models is similar. Figure [I] illustrates the results for model Ml for R m = 10 4 and 
R m = 10 6 and also for model Ml with a reduced magnetic buoyancy effect ( Cp = C v ). The 
developed axisymmetric dynamo regime which is employed at the beginning of evolution 
series shown in Figure [4] is related to the the case R, rn = 10 4 and Cp = 1. This explains why 
the magnetic cycles in model Ml with Cp = C v and R m = 10 6 do not relax to the original 
cycle amplitude. The model with a reduced magnetic buoyancy shows a smaller (by a factor 
2) magnitude of T\\ mode after the relaxation. The mode Tip shows a larger variations of 
amplitude in case of R m = 10 6 than in the case of R m = 10 4 after the relaxation. 

The restoration of the initial of the axisymmetric magnetic field evolution after the 
perturbation does not mean that the non-axisymmetric field completely dissipates. Figures 
g©;a,c ,d) show that a low strength non-axisymmetric field is maintained in the model. 
We find that the strength of the toroidal field T-potentials is reduced from about 100G at 
the maximum to 0.01G after the relaxation. In model M2 the NA magnetic field decays 
slowly after relaxation. Both the Ml and M2 models show a deviation of the parity of the 
AS magnetic field, P from -1 (corresponding to antisymmetric about theequator magnetic 
field) after the perturbation. The NA magnetic field shows a mixed parity solution during 
and after the relaxation phase. The important finding is that the NA perturbation evolves 
through a growing phase for about half an year in model Ml and for about one year 
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in model M2, and it is not just decaying after the initialization. Variations of the axial 
symmetry index, M, and the parity, P, is shown in Fig^a). We see that M first grows and 
then after about 1 year the perturbation is reflected in the parity of the large-scale magnetic 
field. Model Ml shows greater variations of the parities than the model M2. We find that 
in the Ml relaxation of the parity takes more than one cycle. After relaxation the dynamo 
model returns to axial symmetry, M « 10 -6 . Even such a low strength non-axisymmetric 
magnetic field can produce some interesting phenomena which may be related to solar 
observations. 


Figures [5|6^ b) show the longitudinal evolution of a maximum of the large-scale toroidal 
field in the subsurface shear layer. The longitude is computed for the coordinate frame 
rotating with the period of about 25 days. After the perturbation initialization the longitude 
of the m=l mode drifts around the Sun. There is no fluid motion associated with this drift, 
this drift is an analogue of the latitudinal dynamo waves in the slowly-rotating regime. The 


effect was suggested earlier in the study of the linear dynamo regimes by Raedler (1986). 


It was found in observations (e.g., Tuominen et al. 2002b Lindborg et al. 2013) and in the 


direct numerical simulations (Cole et al. 2014). The speed of the drift changes after the 
relaxation. Model Ml shows the almost fixed positions for the longitude of the m=l mode 
during epochs of the maximum of the AS toroidal magnetic field. The latitude-longitude 
position of the maxima of the large-scale toroidal magnetic field strength after relaxation are 
illustrated in Fi^7|(b). Model Ml (without the meridional circulation) shows the periodic 
changes of the longitude by 180° degrees during the magnetic cycle decay. The change of 
the longitude is accompanied by a change of the hemispheric position of the field maximum. 
Thus, the orientation of the global non-axisymmetric field is reversed every cycle during 
the minima of the toroidal magnetic field . This behavior may correspond to the “flip-flop” 
phenomenon of the active longitudes suggested for stellar magnetic cycles (Berdyugina 


2004). 
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Model M2 which includes the meridional circulation has no stable positions of the 
non-axisymmetric magnetic field azimuth. Figures |5|a) and[6](a) clarify the reason for this. 
The effect disappears because the circulation mixes the magnetic field in the subsurface 
shear layer with the magnetic field of the deep interior with a period which approximately 
corresponds to the period of the magnetic cycle. We see that oscillations of the Sjj and 
T\ A harmonics have a 7 t/ 2 phase shift. Thus, for a persistent appearance of the active 
longitude the phases of the S\j_ and Tij harmonics should be consistent. Model M2 shows 
a continuous drift of the longitude of the large-scale toroidal magnetic field strength in 
the course of the magnetic activity evolution. Additional runs which are not illustrated 
here show that the model with the circulation could produce a sort of active longitude 
phenomenon but for another combination of the circulation cells differen from the solar case 
shown in Figur^l^b). More specifically, the strength of the bottom cell should be reduced 
compared to model M2. In this case we find that the active longitudes occupy only the 
Northern hemisphere. 

Figure [8] shows snapshots of the axisymmetric and non-axisymmetric fields after a 
half-year evolution of the initial perturbation. We also show the configuration of the 
external potential magnetic field. This period of time corresponds to a maximum of the 
toroidal magnetic field in the upper part of the convection zone, and the epoch of the polar 
field reversal. The non-axisymmetric part of the field is concentrated to the surface (as 
the initial field). The longitude-latitude diagram shows the distributions of the large-scale 
non-axisymmetric magnetic field. It illustrates how the differential rotation stretches the 
initial magnetic field configuration (cf., Figj2^i). The snapshots in Figures 8 and 9 show the 
large-scale unipolar regions which extend from the equator to the high-latitude regions the 
poles. The increasing of the non-axisymmetric magnetic field in the polar regions results in 
twisted field lines in the polar caps. 
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Figure [TO] illustrates snapshots of the magnetic field configuration during the decaying 
phase of the magnetic cycle in model Ml after 6 years from the initialization. We see that 
the non-axisymmetric toroidal field distributed over the convection zone has maxims at 
the bottom of the convection zones and in the near equatorial region. The strength of the 
non-axisymmetric field is much smaller than the strength of the axisymmetric field. Model 
M2 has long overlaps between the subsequent cycles. Snapshots for this model are shown in 
Figfl0|(bottom) for the growing phase of the cycle. The snapshots show the situation when 
the symmetric with respect to the equator m=l mode dominates at the surface. In the deep 
layers the general distribution of the non-axisymmetric magnetic field is close to model Ml. 

The stationary dynamo evolution begins about 15 years after the initialization of the 
non-axisymmetric magnetic field in model Ml (see Fig[5j) . The relaxation time of model 
M2 is about one cycle. In the stationary stage the non-axisymmetric field is concentrated 
at the top of the convection zone like in the snapshots shown in Figures [8] and [9] (also see 
Figures [5](a) and[6|a)). The antisymmetric m=l mode with mixed parity dominates in both 
models. 


In addition to models Ml and M2 presented in the paper, we calculated the models 
for different initial conditions by changing the spatial distribution of the non-axisymmetric 
perturbation and the initialization time relative to the different epochs of the magnetic 
cycle. In model Ml the effect of perturbation is the greatest when it is initiated at the 
growing phase of the cycle. Also, the impact of the perturbation, and the amplitude 
of the non-axisymmetric field after the relaxation are stronger with the increase of the 


perturbation depth. However, if the perturbation in the form of Eq(28) is located near the 
bottom of the convection zone, it produces only a weak effect on the large-scale distributed 
dynamo. 


To investigate the rotation rate of the non-axisymmetric modes we calculate power 
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spectra of the m = 1 mode azimuth for different radii. Figure 11 shows results for the three 
levels of the solar convection zone. The time series cover the period of about 5 year after 
the initialization of the non-axisymmetric held in the models. At the initialization period 
the equatorial dipole was rotating with a period of 27.26 days, which corresponds to the 
differential rotation period at the latitude of the maxima of the initial perturbation. For 
model Ml it is found that in the subsurface shear layer the equatorial dipole rotates with 
the periods of 25.1-25.5 days. This corresponds to rotation of the subsurface shear layer at 
30° latitude. At the surface the dipole rotates with the period about 25.7 days, and at the 
bottom of the convection zone it rotates with the period of 25.1 days. The origin of these 
rotational periods has to be studied further. It seems that the periods of rotation follow 
the rotation profile in the solar convection zone, e.g., see the typical bow of the magnetic 
held distribution in Fig [Toj'c). In model M2 the meridional circulation mixes all layers of 
the convection zone, producing a unique maximum of the non-axisymmetric magnetic held 
rotating with period of 25.4 days, which corresponds to the differential rotation period at 
latitude of 30° and r = 0.9 R. 


Finally, the animated evolution of the large-scale magnetic held is illustrated by two 
videos: “http://www.youtube.com/watch?v=0i5tslwaxao ’ and 

’"http://www.youtube.com/watch?v=buNK91Sb30A’, showing results for model Ml. 


4. Discussion and Summary 

In the paper we explored the evolution of a non-axisymmetric (NA) magnetic held 
perturbations in the mean-held solar-type dynamo models. The models are kinematic with 
respect to the mean how. The distribution of the mean how is taken from the recent results 
of helioseismology, including the subsurface rotational shear layer and the double-cell 


meridional circulation pattern which was suggested recently by results of Zhao et al. (2013). 
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We studied models with and without meridional circulation. The non-axisymmetric dynamo 
model takes into account the mean turbulent electromotive force in a fairly complete form. 
The mean electromotive force which is employed in the non-axisymmetric part of the model 


is the same as for the axisymmetric (AS) part except that the 5-effect (hi x J term) (Radler 


1969) was omitted in the non-axisymmetric electromotive force. We plan to investigate it 


separately. 


The dynamo models studied in the paper show two distinct evolution phases. The 
transient phase starts after the initialization of the NA perturbation and ends approximately 
after one dynamo cycle period. The results suggest that the non-axisymmetric magnetic 
held can considerably affect the axisymmetric dynamo. This effect depends on the 
amplitude of the perturbation, the depth of perturbation, and the phase of evolution of 
the axisymmetric magnetic held. In the paper we illustrated the effect of a perturbation 
representedby the mix of the odd and even parities of m = 1 mode of magnetic held with 
magnitude of 1G. The depth of the perturbation is r = O.9i? 0 . 

The models show that during the transient phase the NA magnetic held is amplihed, 
to the level of the axysymmetry index M = 0.3, which is 30% of the total magnetic held 
energy. It affects the North-South asymmetry of the magnetic held about equator. The 
growth and decay of the NA magnetic held is accompanied by oscillations with a period 
about of 1.5 year. Rotation of them = lmode during the transient phase shows a continuous 
spectrum of the rotation periods because the evolution of the NA magnetic held is strongly 
coupled with the differential rotation. The maximum in the spectrum of the rotation 
periods is at about of 25.4 days. 


The transient phase of the dynamo evolution follows by a “stationary” phase when 
the NA magnetic held evolves slowly, varying on the time interval which is longer than 
the dynamo period. Solutions for the second phase could be compared with results of the 
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previous NA dynamo models (see, Raedler 1986; Raedler et al. 1990, Moss et al. 1991 


Moss 1999, Elstner & Korhonen 2005). In the “stationary” phase we found only a weak NA 


magnetic field (M ~ 10 6 ) with the dominant m — 1 structure and with the antisymmetric 
relative to the equator magnetic parity. The NA dynamo modes rotates rigidly. 

Model Ml (without the meridional circulation) shows that the rigid rotation of the 
m=l dynamo modes is accompanied by changes of the longitude and the hemispheric 
position of the maximum of the m = 1 mode with a period of one dynamo cycle. This 
regime resembles the so-called “flip-flop” phenomenon which is found in the stellar magnetic 


activity (Jetsu et al. 1991, Korhonen et al. 2002. Berdyugina 2004, Lindborg et al. 2013). It 


was also previously suggested by other NA dynamo models (Tuominen et al. 2002b, Moss 


2004, Elstner & Korhonen 2005 k Berdyugina et al. 2006). The flip-flop phenomenon may 


be related another phenomenon so-called the “active longitudes” (AL) (e.g.. Vitinskii 1966 


Vitinsky et al. 1986). In Introduction it was mentioned that the persistence the AL on the 


Sun on the century time scale remains a highly controversial issue. The origin of flip-flop 
phenomenon in nonlinear non-axisymmetric dynamo models was discussed in details by 


Moss (2004) (see, also, discussion in Berdyugina et al. 2006). In our calculations the 
active longitude is fixed when the T\ \ and S\ t i dynamo-modes are in phase. The flip-flop 
occurs when the orientation of the equatorial dipole changes the sign. If this effect is 
accompanied by the equatorial symmetry variations of the axisymmetric magnetic field then 
the orientation of the large-scale magnetic field changes by 180°. The current explanation 
remains qualitative because the amplitude of the NA magnetic field in the model is rather 
small (the axial symmetry index, M ~ 10 -6 ). It is not clear how such a weak NA magnetic 
field could modulate the sunspot activity. 


One interesting feature of model Ml is that in the stationary evolution phase the 
antisymmetric parity of the m — 1 mode dominates. This is different from previous 
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results of the linear theory by Raedler (1986) and the previous nonlinear models (Raedler 


et al. 

1990; 

Moss et al. 

1991 

Moss 

1999) 


on the mid-latitude zone which is different from the typical poleward concentration of 
the NA magnetic held suggested by results of the above cited papers. The origin of this 
disagreement is unclear. We have to stress that in our models we employ a fairly complete 
information about the radial profiles of the a-effect and the differential rotation suggested 
by the modern results of the helioseismology, and results of recent theoretical works. Also 
we include the magnetic helicity conservation which was not taking into account before in 
the non-axisymmetric dynamo models. The dynamo parameters such as the a-effect and 
the turbulent diffusion coefficients were tuned to reproduce the 22 year dynamo cycle by the 
axisymmetric dynamo. Considering the ratio between the a and ^-effects in our models, 

1 1Q I // j n 

we find -- ~ 10 , where a$ is the magnitude of the a-effect. This is about by 2 orders 

ao 


of magnitude than the quantity employed by Raedler et al.| ( |1990 ) and Moss (1999). The 
finding similar to ours, i.e., a concentrated to the equator antisymmetric m = 1 magnetic 


field mode was reported by Nelson et al. (2013) for the global numerical simulation of the 


convective dynamo on the solar-type star rotating 3 times faster the modern Sun. 

Model M2 with the meridional circulation illustrates another feature predicted by the 
linear analysis of the kinematic NA dynamo models. In this model the longitudinal position 
of the m = 1 mode in the stationary phase of evolution drifts around the Sun with a period 
of about one dynamo cycle. This can be interpreted as an azimuthal dynamo wave. The 
effect was predicted by the mean-held theory (see, e.g. Krause & Radle r|[l9 80), it was found 
in the mean-held NA dynamo models (Raedler] 1986; Raedler et al. 1990; Moss et al. 1991), 


and in the direct numerical simulations (Cole et al. 2014). It was also found in observations 


of magnetic activity on fast-rotating late-type stars (e.g., Tuominen et al. 2002b; Lindborg 


et al. 2013) 
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5. Conclusions 


We considered non-axisymmetric mean-field dynamo models including the non-linear 
magnetic helicity and magnetic buoyancy effects, which were not studied before. The 


study confirms the previous findings of Moss (1999) that the non-axisymmetric dynamo 
component is rather weak if we start from a weak initial (seed) non-axisymmetric field. We 
notice that our models can be characterized as weakly non-linear because the parameter 
of the o-effect in the model is only 30% above the dynamo instability threshold. Also 
the magnetic helicity conservation and magnetic buoyancy prevent the generation of 
magnetic field of the super-equipartition strength (|Brandenburg &; Kapyla 2007, Hubbard 


& Brandenburg 2012). Thus, the low-strength non-axisymmetric magnetic field generated 


from a weak seed field can be explained by the linear stability of the non-axisymmetric 
field, and by the weak non-linearity of the dynamo system. However, finite-amplitude 
non-axisymmetric perturbations, which can be developed in the complex dynamical system 
may have significant effects on the dynamo process. 


The modeling results show that the magnetic helicity conservation (also known as 
dynamical quenching of the a-effect) is an important factor to preserve the non-axisymmetric 
field from a complete decay. It is found that for the magnetic Reynolds number R m = 10 6 
the coupling of the non-axisymmetric magnetic field with the axisymmetric dynamo process 
is stronger than in case of R m = 10 4 . The strong coupling results in a synchronization 
between oscillations of the non-axisymmetric and axisymmetric magnetic fields. Also, our 
models include nonlinear effects of magnetic buoyancy. It is found that if the magnetic 
buoyancy effect is switched off then the strength of the non-axisymmetric field after the 
relaxation is decreased by a factor of two. 


The paper illustrates our initial results of the nonlinear non-axisymmetric mean-field 
dynamo model. The axisymmetric part of the model is based on our previous results. The 

















28 


non-axisymmetric perturbations are assumed to be located in the near-surface rotational 
shear layer. For the first time we demonstrate that nonlinear coupling between the 
asymmetric and the non-axisymmetric fields can impact the generation of the axisymmetric 
held in the case of finite-amplitude perturbations. The effect depends strongly on the 
dynamo mechanisms involved in the problem, the spatial distribution of perturbation, the 
phase of the dynamo cycle at the time of initialization of the perturbation, and on how well 
the magnetic helicity is conserved in the system. These factors determine the subsequent 
evolution of the dynamo system including the dynamo cycle, evolution and rotation of the 
non-axisymmetric modes of large-scale magnetic fields. 

In summary: 

• The differential rotation and interaction with the axisymmetric magnetic fields results 
to amplification of the non-axisymmetric perturbation during the transient phase of 
evolution, producing the magnetic held configuration which deviates considerably 
from the axial symmetry. 

• In the solar dynamo models non-axisymmetric magnetic held perturbations developed 
in the near-surface rotational shear layer affect the strength and period of the dynamo 
cycles. 

• Without the meridional circulation the non-axisymmetric dynamo-mode shows a 
“hip-hop” phenomenon and also has a hxed longitudinal position resembling the active 
longitude phenomenon. However, the solar-type (double-cell) meridional circulation 
destroys these effects. Instead, the non-axisymmetric held represent a travelling in 
azimutal direction dynamo wave. 
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Fig. 4.— Evolution of the simulated sunspot number (a) and T\_\ harmonic (b) in model Ml 
for the different values of the parameter R m and Cp: R, rn = 10 4 , Cp = 1 (red), R rn = 10 4 , 
Cp = C v (blue), R m = 10 6 , Cp = 1 (black). 
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Fig. 5.— Non-axisymmetric modes of model Ml: a) the time-radius evolution of the dynamo 
modes T M (background image) and Si t i modes (contours are in the same range of values as 
the color scale); b) evolution of longitude of the maximum of the large-scale toroidal field, c) 
the mean strength of the near-surface axisymmetric and non-axisymmetric toroidal magnetic 
field and the strength of the m=l T-potentials, d) evolution of parities, P and P, and the 
index of the axisymmetry M. 
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Fig. 6.— The same as Figure [5] for model M2 
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Fig. 7.— a) Evolution of parity of total magnetic field P versus the axisymmetry index M 
for model Ml (blue) and M2(red); b) Position of maxima of the total toroidal magnetic held 
strength of the subsurface shear layer (r=0.92R) in models Ml and M2 after relaxation. 
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Fig. 8.— Model Ml a half year after the initialization of the non-axisymmetric perturba¬ 
tion, a) distribution of the axisymmetric toroidal magnetic field (background image) and 
poloidal field lines in the meridional cross-section, b) azimuth-latitude distribution of the 
non-axisymmetric radial magnetic field (background image) and toroidal magnetic field (con¬ 
tours are in the same range as the color-scale) on the surface; c) the same as in panel (a) 
for the non-axisymmetric components of magnetic field; d) the magnetic field lines of the 
non-axisymmetric magnetic field above the surface. 
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Fig. 9.— The same as Figure [8] for model M2 
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Fig. 10.— The same as in Fig. 7 after 6 years of the evolution after the initialization of the 
non-axisymmetric perturbation. 























































RADIUS 


37 



Fig. 11.— The power spectra for rotation periods of azimuth of non-axisymmetric mode 
,S'i ) i in the solar convection zone for: a) model Ml; b) model M2. 
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